From blockade to transparency: controllable photon transmission through a circuit QED system 
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A strong photon-photon nonhnear interaction is a necessary condition for photon blockade. Moreover, this 
nonlinearity can also result a bistable behavior in the cavity field. We analyze the relation between detecting 
field and photon blockade in a superconducting circuit QED system, and show that the photon blockade cannot 
occur when the detecting field is in the bistable regime. We further demonstrate that the photon transmission 
through such system can be controlled (from photon blockade to transparency) by the detecting field. Numerical 
simulations show that our proposal is experimentally realizable with cuiTent technology. 

PACS numbers: 42.50.Gy, 42.65.Pc, 42.50.Ar, 85.25.-j 

Introduction. — Superconducting circuit quantum electrody- 
namics (QED) allows studying the interaction between artifi- 
cial superconducting atoms and quantized microwave fields 
(see, e.g., reviews fl, 7] and references therein). The cou- 
pling strengths have been explored from the strong regime to 
the ultrastrong one |l3|-|6J]. It is well known that the equally- 
spaced energy structure of the quantized microwave field can 
be changed to an anharmonic one, including dressed states 171], 
formed by a superconducting charge qubit and classical mi- 
crowave field, which have been observed experimentally fS]. 

The nonlinear energy splitting in circuit QED has been 
shown by experimentally analyzing the Rabi frequencies for 
different numbers of microwave quanta inside the cavity [9]. 
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Moreover, the experimental spectrum for the square root of 
the photon-number nonlineaiity was reported llioll . Further- 
more, the nonlinear response of the vacuum Rabi resonance 
was demonstrated [11]. CuiTent experimental data indicate 
that the nonlinearity of the microwave photons can be many 
orders of magnitude larger than that of macroscopic media. 
All of these experiments I.8. -.1 1.1 lay a solid foundation for the 
development of microwave nonlinear interactions on a chip, 
and open the door to further study microwave nonlinear quan- 
tum optics at the level of single artificial atoms and single mi- 
crowave photons. For example, photon blockade [12], which 
is the analog of the Coulomb blockade for electrons in single- 
electron devices, has recently been observed in: an optical 
cavity with one trapped atom lll3ll . quantum dots coupled to 
the photonic crystal cavity 1 14], and also superconducting cir- 
cuit QED with resonant [15] and dispersive 1.16.1 atom-field 
interactions. 

Photon blockade originates from the photon-photon inter- 
action, which prevents subsequent photons from resonantly 
entering a cavity. Therefore, photon blockade can be used as 
a single-photon source or single-photon transistor. The mi- 
crowave nonlinear interaction might also be used to improve 
superconducting qubit readout ]17, 18]. One of the basic 
conditions for photon blockade is that the decay rate of the 
cavity field should be less than the photon-photon interacting 



FIG. 1: (Color online) (a) Eigenvalues Euin^ (green curves) of the 
free Hamiltonian hK(n — ric) versus ric for the photon numbers n = 
and n — 1 inside the cavity, in energy units of tin. Blue curves are 
the eigenvalues near ric ~ 0.5, when the external field is included. 
The degeneracy of the two eigenvalues in the free Hamiltonian is 
lifted at the point ric = 0-5, with distance h\i}\ by the external field, 
(b) Several eigenvalues Ek{nc) (up to four) of the Hamiltonian in 
Eq. (O vary with ric for the ratio K/\fl\ = 10. (c) The mean photon 
number in the ground state varies with ric for n/\n\ = 1 (green dash- 
dotted curve), 10 (red dashed curve), and 100 (blue solid curve). This 
is the photon analog of the Coulomb staircase. 



strength. We know that the decay rates of the photon number 
states are very different from those of coherent states 1^191 12011 . 
thus the blockaded photons might not be observed when the 
photon number inside the cavity is extremely large. In addi- 
tion, bistability is one of basic properties of nonlinear sy stems, 
but there is lack of studies of bistability on the photon block- 
ade. 

Here we focus on the Jaynes-Cummings model of circuit 
QED in the strong-coupling regime, study the relation be- 
tween bistability and photon blockade, and analyze the effect 
of the detecting field on the photon blockade by studying the 
second-order degrees of coherence for the cavity field, explore 
the possibility of unblocking the photon of the weak detecting 
(or signal) field by using electromagnetically-induced trans- 
parency ]122(]. 

Model. — We start from the standard Jaynes-Cummings 
model for circuit QED with a transmission line resonator cou- 
pled to a superconducting qubit with coupling strength g. The 



TABLE I: Equivalence between photon and Coulomb blockades 
(or similarity between photon blockade and superconducting charge 
qubits). 



Quanta 


Photons 


Electrons 


Characteristic energy 


Ken- energy flK, 


Charging energy E^ 


Energy offset control 


Driving field frequency uJd 


Gate voltage Vg 


Detecting method 


Driving field £;{0 


Bias voltage Vt 


Measurement of output 


Photon number (71) 


Electron current / 


Tunnehng energy 


Coupling energy h\il,\ 


Josephson energy Ej 



frequencies of the cavity field and the qubit are luq and tUq, 
respectively. If the cavity is driven by classical fields, then 
with a canonical transformation Bill , the driven Hamiltonian 
is transformed to the dressed one 
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FIG. 2: (Color online) Logarithm log(|As|) of the steady-state so- 
lution \As\ versus the strength of the driving field \Q\, for different 
detunings in (a) and (b). Here, A^up — g'^/A'^'y = K/7 is the up- 
per bound photon number for the photon blockade, A/2tt — 1 GHz, 
7/27r = 0.1 MHz, and g/2-K = 200 MHz. For above given pa- 
rameters, \As\ has two stable states when 0.0091 MHz < Ad < 
36.95 MHz. 



H = hujQa^a H — [ujq 



E{N)] a, + h [E{t)a) + E*{t)a] ,(1) 



with the creation a^ (annihilation a) operator of the cav- 
ity field and E{N) = ^A^ + Ag'^N. Here the operator 
N = a^a + {az + l)/2 describes the total number of exci- 
tations for the photon and the qubit, and A — ujq — ujg de- 
scribes the detuning between the cavity field and the qubit. 
The third term in the right of Eq. ([T), with the parameter E{t), 
describes the interaction between the external fields and the 
cavity mode. In the derivation of Eq. ([T]i, the terms proportion 
to (9(iV~^/^) are neglected. The excited and ground states 
with sign ± for az are considered as the attractive and repul- 
sive photon-photon interactions. Below we assume that the 
qubit is always in its excited state. We also assume that the 
decay of the qubit is negligibly small in the bad-cavity limit. 
If a monochromatic driving field (or saying the detecting 
field for photon blockade), with frequency ujd, is coupled to 
the cavity via E{t) = J7e^'"''* with coupling constant il, 
then the Hamiltonian in Eq. ([TJ can be used to describe pho- 
ton blockade [,12,] when the amplitude of il for the driving 
field is much smaller than the photon-photon coupling con- 
stant K given below. Moreover the decay rate 7 of the cavity 
field should also be much smaller than k to experimentally 
observe photon blockade. In the dispersive regime, the single 
two-level atom-induced photon blockade can be understood 
by expanding the Hamiltonian in Eq. ([T]i up to third order in 
parameter g/A as 



Hcs — hnin — Tie) + h{Vta* + 0*a) 



(2) 



with n — at a in the rotating reference frame with frequency 
LOd for the driving field. Here, k = g'^/A^ denotes the photon- 
photon interaction strength, Uc = (k — Ad) /2k is the rescaled 
detuning A^ — ujq — uJd between the driving and the cav- 
ity fields. In the derivation of Eq. (|2]l, several constant terms 
and also the qubit state-dependent cavity-frequency shift have 
been neglected with the assumption that the qubit is in its ex- 
cited state. Equation (|2]i shows that Hk and tud (or Uc) are 
equivalent to the charging energy and gate voltage in single- 
electron devices for Coulomb blockade. Equation (|2) also 



shows that the parameters k, tUd (or jic), and \il\ play simi- 
lar functions as the charging energy, the gate voltage, and the 
Josephson energy of the charge qubit (see, e.g., Ref. 112311 '). For 
resonant driving, A^ = and ric — 0.5, the photon states |0) 
and |1) are degenerate for the free Hamiltonian hK{n — 0.5)^. 
As shown in Fig. \lia), this degeneracy can be lifted by the 
coupling strength \n\. Figures [Tfb,c) show the variations of 
the eigenenergies and the mean photon number (n) of the 
ground state for the Hamiltonian in Eq. (12} as functions of 
the parameter ric- Variations for both the eigenenergy and the 
mean photon number are the same as those of the eigenenergy 
and mean charge number of charge qubits or single-electron 
devices for Coulomb blockade. The staircase shape for the 
mean photon number was experimentally demonstrated in cir- 
cuit QED system fll6ll . Moreover, Fig. [Ha) also shows that a 
large ratio K/\il\ corresponds to a sharper step, i.e., a weak 
driving field and strong photon-photon interaction are more 
useful for detecting photon blockade. Photon and Coulomb 
blockades are compared in Table I. 

Bistability and blockade. — As shown in Fig. [Uc), the basic 
feature of the staircase shape for the photon blockade gradu- 
ally disappears when the strength of driving field approaches 
or is larger than the photon-photon coupling strength. It 
is known that the driven nonlinear photonic systems exhibit 
bistability in a regime of parameters. To study the relation 
between bistability and photon blockade, let us now write the 
following Heisenberg-Langevin equation 



da 



iu}Q + "f ~ i 



g^cTz 



EiN) 



ifle^"^"* 



27ai„(i)(3) 



for the driven Hamiltonian in Eq. ([T) by the monochromatic 
classical field E{t) = J7e~*'^''*. The parameter 7 is the de- 
cay rate of the cavity field. The quantum fluctuations ai„ 
due to the vacuum field are assumed to have a zero mean 
value (ain(i)) ~ and Markov coiTelation {ainit)al^it')) — 
d{t - t'). If we denote a{t) = A(t)e''-'^''\ then the steady- 
state solution As = {A)s for the cavity field becomes 



As 



-in/{j+i[Ad-g^<Jz/EiN)]}, 



(4) 



by using Eq. Q and the mean-field approximation, e.g., 
{a^a)s = {a^)s{a)s- Here E{N) is the expression of E{N) 
with replacing operator a^a in A^ by steady state value |^sP- 
Equation (01) shows that the cavity field amplitude |Asp 
has two stable states when A^ satisfies the condition A_ < 

Ad<A+withA± = [gyE(N)]-{2g*\A,\^T[i9'\As\*- 
^^ E^{N)]^/^} / E^{N)]. Equation © also clearly shows 
when the qubit is decoupled from the cavity field (i.e., g — 
0) or when the driving field makes \As\'^ extremely large, 
such that gaz/E{N) « 0, then the response of the sys- 
tem is the same as that of the harmonic oscillator How- 
ever when g ^ 0, the resonant peak of \As\^ will move to 
Wfj = cjq — 9<yz/ E{N). That is, when the cavity contains 
\As\'^ photons, the frequency LOd of the driving field should 
be increased an amount —gaz/E{N) to overcome the photon 
blockade [12]. The upper bound photon number inside the 
cavity for the photon blockade is A^up ^ gtJl^'^l) — '*/7 ^^ 
the dispersive regime as discussed in Ref. MM- 

In Fig.|2l the steady state \As \ versus the input \il\is plotted 
for several different detunings A^. Figure|2]clearly shows that 
the bistability disappears for either A^ > A_)- or A^ < A_. 
Figure|2]also shows that most values for \As\ are smaller than 



the upper bound value 



for some values of A^ in the 



bistable regime, e.g., values in the region near Ad/2TT — 37 
MHz, as shown in Fig. I2b). However, as shown in Fig. |2] 
(a) for other parameter regimes of Ad, we find that the up- 
per bound value ^ 



up can be smaller than some values of 
\As\ corresponding to the lower branch of the bistable curve. 
Therefore, Fig. |2]tells us that A'^^ip is a necessary but not suffi- 
cient condition for photon blockade. Because one input corre- 
sponds to two stable outputs in the hysteresis region, thus the 
photon blockade is not well defined in such region. 

Dependence of photon blockade on the detuning Ad and 
coupling strength Vt. — To show the effect of the driving field 
on the photon blockade, let us now study the statistical proper- 
ties of the cavity field when the input Oin of the vacuum fields 
in Eq. ^ is considered. We assume that the vacuum fields Oin 
in Eq. (|3]l result in a small fluctuation Aj [t] of the cavity field 
near its stable steady-state Ag by writing the cavity operator 
as A{t) = As+ A fit). Besides the steady-state solution as in 
Eq. (|4]l, with the input ain(i) — Ain(i)e^'"''', we can obtain 
an equation of motion for the fluctuation operator Af [t) as 



dAf{t) 



dt 



iAd{N)+-i 



Af (t) -ia{N)A} {ty^2jAin{t) 

_ _ _ 15) 

with AdiN) = Ad - 6uj{N), a{N) = 2g'^A^^z/E^{N), 
and Suj{N) = g^az[E^(N) - 2g^\As\^]/E^(N). Where 
the terms with higher orders of Af{t) and A^{t), e.g., 
the term Af{t)Af{t), have been neglected when Eq. (|5]l 
is derived. By applying the Fourier transform Af{t) = 
f^ e~'^'^^Af{uj)duj/^/2TT and also using the conjugate of 

Eq. (|5]l with At{uj) — [Af{~-u})]^, we can easily obtain the 
solution of the fluctuation as 



Afico) 



V27 
diu}) 



(Ad - w + ij^ Ai„(a;) + a(7V)4n(w) 
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FIG. 3: (Color online) p'^' (0) versus the driving field strength \Q,\ in 
(a) and g^^^ (r) versus the delay time r in (b). Here the other param- 
eters are taken as Ad = 38.5 MHz, A/27r = 1 GHz, g/27r = 200 
MHz, and 7/27r = 0.1 MHz. g*^'(0) versus the detuning Ad 
in (c) for g/27r = 100 MHz and (d) for g/2TT = 200 MHz; 
with A/27r = 1 GHz, l!:!l/27r = 0.01 MHz, and 7/27r = 0.1 
MHz. For parameters given in (c) and (d), the condition of stable 
states for the driving field is 0.018 MHz < Ad < 9.627 MHz and 
0.0091 MHz < Ad < 36.95 MHz, respectively. 



with the denominator factor d(w) = (ia;+7)^+A^— |a(iV)p. 
The statistical properties of the cavity field can be described 
via the second-order degrees of coherence g'^'^^r). Because 
we assume that the vacuum input is Gaussian and satisfies the 
Markov coiTelation, thus g^'^^r) can be obtained by just cal- 
culating the coiTelation function with two operators using the 
Wick's theorem, that is 



9^'Hr) 



2MAl{A\{t)A){t')) + \A,\\A\{r)Af(t'))] 

{\AsV + {nf{tW 
\{A\{t')A].{t))\' + \{A]{t')Af{t))\' 



{\As 



{nf{tW 



1 (7) 



(6) 



with f' = f + T and (n/(t)) = (Aj.(i)A/(i)). It is easy to 
obtain g^'^^r) straightforwardly by calculating all coiTelation 
functions in Eq. ^ using Eq. (|6]l and [Af{—uj)X^. The second- 
order degrees of coherence .g^^'(O) and ^'■^•'(t) are plotted in 
Figs.|3ja,b) using experimentally-accessible parameters. Fig- 
ures |3ja,b) show that the cavity field tends to the classical 
behavior when increasing the strength Vl. In this case, the 
photons might not be blockaded and can transparently pass 
through the circuit QED system. To explore the effect of the 
frequency of the driving field, g'^'^^Q) is plotted as a func- 
tion of the detuning A^ in Figs. |3|c,d) for different coupling 
strengths between the qubit and the cavity field with other pa- 
rameters given in the caption of Fig. |3] Figures |3lc,d) clearly 
demonstrate that the nonclassical behavior of the cavity field 
is out of the bistable regime for detuning A^. For example. 
Fig. Oc) shows g^'^\Q) < 1 when Ad/27r > 9.85 MHz, 
which is larger than the upper bound value 9.627 Mhz of A^ 
for the bistabilility, and thus the photon blockade cannot occur 
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FIG. 4: (Color online) The real Ar(lj') and imaginary Ai{uj') parts 
as a function of cj' = iUd ^ co are plotted in (a) and (b), respectively, 
for two different values of A^. 



in the bistable regime. 

Photon blockade and transparency. — To further discuss 
properties of the Ught field transmission when the photons 
are blockaded, we now study the response of the circuit QED 
system to the vacuum input field, which can be considered 
as a weak detecting field. According to the input-output the- 
ory ll28ll . the output field with the Fourier components can be 
written as 



aout(w) = Ain{uj') + ^/¥fAf{uJ') 



27 



2-fA,+in' 

5{uj), 

(8) 

with uj' = CO + ojd and Af{uj) given in Eq. (|6]l. The physical 
meaning becomes clear if the vacuum input ain{t) is assumed 
as a single-mode field ain(t) — ^exp(— iwi) with the real 
parameter f <C |^|; that is, this weak detecting field does not 
change the statistical properties of the cavity field. In this case, 
uj' in Eq. (|8]l is changed to ui' = ujd — ^, and the terms Ajn (w') 
and 6{uj') in Eq. ([8]l denote the response of the system to the 
input vacuum field and the classical driving field, respectively. 
However, the term with Aj^^{uj') exhibits four-wave mixing 
with frequency uj — 2ujd, which will be studied elsewhere. 

The coefficient of Ain{uj') for the output in Eq. (O corre- 
sponds to that of Ain{uj') in the expression Af{uj') of Eq. (|6]l 
plus one. Thus in Figs. |4ta,b), the real A-r,{lj') and imagi- 
nary Ai{lu') parts for the normalized coefficient of Ain(i^)/C 
in Eq. (|6]l for Ain{t) = ^exp[— i(u; — U!d)t] are plotted by 
using the same parameters as in Fig. |3jc), with A^ « 9.96 
MHz and A^ w 9.74 MHz. These two values correspond 
to the minimum (photon blockade) and the maximum (clas- 
sical case) of ^'^'(O) in Fig.[3jc). As expected, we find that 
the response of the circuit QED system to the input field (or, 
say, weak-detecting field) has a Lorentzian shape, which is 
the same as the decay spectrum of the single photon, when 
the photon is blockaded. However, the weak detecting field 
shows transparency windows to the circuit QED system when 
the driving field is changed such that the cavity field is in the 
classical regime (or say the photon is not blockaded). Thus, 
we can control the photon (from blockade to transparency) by 
changing the applied classical field. 

Conclusions. — In summary, we have studied the tunable 
transmission from the photon blockade to the photon trans- 
parency in superconducting circuit QED systems when the in- 
teraction between the qubit and the cavity field is in the dis- 



persive regime. We analyze the effect of the driving field on 
the photon blockade and show the relation between the optical 
bistability and the photon blockade. We find that the photon 
blockade can be controlled by a classical driving field, that is, 
the photon blockade strongly depends on the properties of the 
driving field. We also find that the circuit QED system can be 
used to generate a four-wave mixing signal. All parameters in 
our numerical simulation are taken from experimentally avail- 
able data. Therefore, our study should be experimentally real- 
izable with current technology. 
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